Theory of monolayers with boundaries: Exact results and Perturbative analysis 
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Domains and bubbles in tilted phases of Langmuir monolayers contain a class of textures known 
as boojums. The boundaries of such domains and bubbles may display either cusplike features 
or indentations. We derive analytic expressions for the textures within domains and surrounding 
bubbles, and for the shapes of the boundaries of these regions. The derivation is perturbative in 
the deviation of the bounding curve from a circle. This method is not expected to be accurate 
when the boundary suffers large distortions, but it does provide important clues with regard to the 
influence of various energetic terms on the order-parameter texture and the shape of the domain 
or bubble bounding curve. We also look into the effects of thermal fluctuations, which include a 
sample-size-dependent effective line tension. 
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I. INTRODUCTION 

Monolayers of insoluble surfactant molecules confined 
to the air/water interface possess complex phase struc- 
tures In the "tilted" phases, the long axes of the sur- 
factant molecules in the monolayer are uniformly tilted 
with respect to the normal and the molecular tilt azimuth 
organizes spontaneously on macroscopic length scales. 
The structures adopted by the molecular tilt azimuth 
are referred to as textures. There is no long-range order 
of the tilt azimuth in the liquid expanded (LE) and the 
gaseous (G) phases. Tilted phases can coexist with the 
LE and G (isotropic) phases and form micron-sized do- 
mains. Alternatively, bubbles of an isotropic phase may 
appear against a background of a tilted phase. Nontrivial 
textures in the domains, and around the bubbles, have 
been observed in the L2/LE and L2/G coexistence re- 
gion, where the L 2 phase is one of the tilted phases. Boo- 
jum textures, similar to those seen in superfiuid 3 He J2] 
and smectic-/ (Sm-/) droplets in liquid-crystal films jl, 
have been observed in the interior of L2 domains . An 
"inverse boojum," which is the texture around the bub- 
ble analogous to the boojum in the case of the domain, 
has been identified ^ . The domains and bubbles associ- 
ated with boojums are not circular. Among the nontriv- 
ial domain shapes seen are protrusions on both bubbles 
and domains, at times sharp enough to be characterized 
as "cusps" iQpJ and indentations in domain boundaries 
which impart a cardioid appearance to the domain [BLQ . 
Such domains and bubbles with unusual textures and 
shapes ought to be observable in other tilted phases as 
well. 

The above textures can be understood in terms of con- 
tinuum elastic theories of smectic liquid crystals ^ . The 
bulk energy is controlled by elastic moduli that quantify 
the energy cost of bend and splay distortions. There are 
also contributions to the boundary energy, known as the 
line tension, that depend on the relative angle between 
the boundary normal and the tilt azimuth. In equilib- 
rium, the texture in a domain or surrounding a bubble, 



and the shape of the boundary between condensed and 
expanded regions, adjust so as to minimize the total en- 
ergy of the monolayer. Domains with nontrivial textures 
and shapes represent the compromise arising from the 
competition between the bulk energy and the line ten- 
sion. 

Simultaneous determination of the texture and the 
boundary of the domain poses a calculational challenge. 
Earlier studies include the exact result discovered by 
Rudnick and Bruinsma for a domain with isotropic elas- 
tic energy and only the first anisotropic contribution 
in the Fourier expansion of the line tension a{4>) ||, 
and the perturbation about the exact result in terms 
of coefficient of the second anisotropic term in the ex- 
pansion |)J. Galatola and Fournier have approached 
the problem of domains with elasticity and line-tension 
anisotropies by searching numerically for the equilibrium 
domain shapes a nd p ositions of domains in a fixed tex- 
ture background ||lof . Riviere and Meunier [|J have at- 
tributed their experimental findings on domain shapes 
and textures to elasticity anisotropy in the same manner 
as in Ref. JlO| |. In the work of Fang et al. 0], nontriv- 
ial boundary shapes for both domains and bubbles, as 
well as the "inverse boojum" textures in the L 2 phase 
surrounding the bubbles, have been reported. A brief 
account of the theoretical understanding of the bubbles 
has also been presented in Ref. ||. 

In this work, we extend the effort of Rudnick and 
Bruinsma to analyze the problem of domains with 
anisotropic clastic energy in addition to the line-tension 
anisotropy. We also generalize the approach to the prob- 
lem of bubbles and provide a detailed derivation of the 
results that have been published in Ref. j|] . Careful anal- 
ysis reveals that although protrusions can be expected to 
form on the boundary of a domain of the L 2 phase, a 
"cusp" in the form of a discontinuity in the slope of the 
bounding curve surrounding the domain will not appear 
in the parameter regime that is appropriate to the anal- 
ysis that has been carried out [p|pX|j5|] . The conclusion in 
Ref. P that a cusp exists is due to an approximation 



1 



that affects the qualitative results of the analysis. The 
fact that the cusp does not exist and the condition for the 
existence of cusps on the boundary were first pointed out 
by Galatola and Fournier |Io| . A formal derivation of the 
conditions for the appearance of a cusp will be provided 
in this paper. Perturbative results, which yield the effects 
of small anisotropics on the textures and boundaries, are 
obtained. The reliability of the perturbative approach 
when the boundary is significantly different from a circle 
is not obvious. Nevertheless, one is provided with useful 
insights with regard to the influence of various contri- 
butions to the energy of the Langmuir monolayer. In 
addition, we examine the effect of thermal fluctuations. 
We are led to a renormalized line tension that depends 
on the radius of the boundary Jll|] . 

We have also implemented a numerical program using 
finite element methods for evaluation of the equilibrium 
texture and boundary simultaneously. With the use of 
this approach, we are able to explore regions of the pa- 
rameter space that are not accessible to the perturbative 
technique. A brief report on the numerical work has al- 
ready appeared p3| . A full description of this method 
and a systematic review of the results of its implementa- 
tion are deferred to a future article. 

The organization of this paper is as follows. In Sec. 
||, we describe the approach in general. In Sec. |l|, we 
summarize the exact analytic results. Section ^ displays 
the perturbative analysis of the relation between the tex- 
ture and the boundary. Sections |y| and |Vl] describe the 
analysis for the cases of domains and bubbles that results 



from perturbing about the exact solutions. In Sec. VII 



we analyze the effect of thermal fluctuations. Concluding 



remarks are contained in Sec. VIII 



II. THE APPROACH 

We describe the monolayer by an ordered phase 
with XY"-like order parameter c(x,y) = xcosQ(x,y) + 
y sin 9(cc, y), a two-dimensional unit vector indicating the 
direction of the projection onto the substrate of the tilted 
hydrophobic tail of the surfactant forming the Lang- 
muir monolayer. The quantity, Q(x,y), is the angle that 
c(x, y) makes with the x axis. When a region Q contains 
an ordered phase which is invariant under in-plane re- 
flection, the free energy of the system takes the general 
form § 

H[Q{x,y)}= { H b dA+ I a[3 - Q{x,y)]ds, (2.1) 
Jq Jr 



where 



H b = ^|V • c(x, y)\ 2 + ^ |V x c(x, y)| 2 , (2.2) 

(2.3) 



<r((f>) = (Jq + a n cosnq 





FIG. 1. The geometry of the calculations for (a) domains 
and (b) bubbles in plane-polar coordinates where the bound- 
ary r is parametrized by p((p). The gray area is the bulk 
designated by Q. h and t are the outward normal and the 
tangent, respectively. is the angle between the c director 
and the x axis and -d is the angle between the outward normal 
of the boundary and the x axis. 

elastic moduli, and $ is the angle between the outward 
normal of the boundary and the x axis. The quantity 
(To > is the isotropic line tension. The first integral is 
over the area, fl, of the system, while the second is over 
the boundary, T, as indicated. The setup of the problem 
in plane-polar coordinates is shown in Fig. [lj. 

Minimization of the energy leads to equations for 
Q(x,y) and the bounding curve T. Q(x,y) satisfies 



Qyy) COS 26 



29^ sin 29 



b[(G xx 

(-el + Ql) sin 29 + 29^9,, cos 29] = (2.4) 



in Q and 



Here, K s and Kb are, respectively, the splay and bend 



K9„[l-6cos2(tf-9)] + 
K be t sin 2(tf - 9) - <r'(d - 9)= 



(2.5) 



2 



along r, where 9„ = n ■ V9, 9f — i ■ V0, h and i are, 
respectively, the outward normal and tangential vectors, 
and 



K s + K b 



K s — K b 
K« + Kt, 



(2.6) 
(2.7) 



The primes attached to functions denote derivatives, e.g., 
<j'(<fi) = da((f>)/d(f>. The extremum equation for the 
bounding curve T, implicitly in terms of O n , 9t, and 

d-d/ds, is 



H b -a'(#-&)e n -a"(ti-@)G t 

+ [a(6 - 9) + a"ifi - 6)] ^ + A= 

ds 



0, 



(2.8) 



where ds is the length element of Y traversing in the pos- 
itive direction of Q and A is a Lagrange multiplier that 
enforces the conditi on o f cons tant e nclo sed area. The set 
of equations, Eqs. (2.4), (2.5), and (2^8), are highly non- 
linear. It appears, in general, impossible to find general 
analytical solutions to this set of equations. However, 
there arc full analytical solutions for special cases. 




FIG. 2. The c director distribution and the BAM re- 
flectance in a domain computed for k = 1, Ro = 20, and 
a v — 1.6, where p = 1 in (a), p — 2 in (b), p — 3 in (c), and 
p = 4 in (d). 



III. EXACT SOLUTIONS 

We start with the assumption of a circular boundary. 
We restrict our considerations to isotropic elastic moduli, 
i.e., 6 = 0. Additionally, we assume that the a niso tropic 
line tension, as given by the expansion in Eq. (2.3), con- 
tains only one term, in that a n — for all n ^ p. We 
will take a p > 0. This is because the texture with a p < 
can be trivially obtained by rotating all c{x, y) simulta- 
neously by (2m + l)7r/p, where m is an integer S [0,p— 1], 
due to the symmetry in the line tension. In this special 
case, Eq. (|2.4|) reduces to Laplace's equation 



V 2 6 



0. 



(3.1) 



and Eq. (2.5) in the plane-polar coordinate system be- 
comes 



k<3 p — a'(ip - 
k@ p + a'(n + tp ■ 



9) 
9) 



0. 
0. 



(3.2) 
(3.3) 



wher e Eq. ( |3.2| ) applies to the case of a domain while Eq. 
( |3.3| ) is appropriate to the case of a bubble. In two dimen- 
sions, the solution to Laplace's equation can be written 
in general as 
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9(fc,^) = -[/(e fe+ ^)-/(e fe -^)], 



(3.4) 



with f(z) an analytic function of z = e k+i(p in the region 
of interest, £1 for our case. In the case of a circular do- 
main of radius Rq centered at the origin, it is shown in 
Appendix ^| that 



fo(z) 



- ln(l — a p z p ), 
P 



= -e+^/l + e 2 , 



(3.5) 
(3.6) 
(3.7) 



satisfy Eq. (3J2). We have defined here a dimensionless 
parameter e = n/(pa p Ro). Figure [2] illustrates such solu- 
tions for p = 1, 2, 3, and 4. Also displayed in the figure 
on the background of each plot is a simulation of the 
image that would be obtained by Brewster angle micro- 
copy (BAM). The BAM reflectance depends on the exact 
experimental setup and the properties of the monolayer. 
A detailed discussion on the computation of the BAM 
reflectance can be found in Ref. p| . In the case of all 
simulated images presented in Fig. ^ and elsewhere in 
this paper, the Brewster angle is taken to be that of wa- 
ter Ob = 53.12°, the angle of the analyzer a is equal 
to 90°, the thickness of the monolayer is assumed to be 
d = 0.3 ?im, the tilt 'F is 30°, the dielectric constants 
of the monolayer are e± — 2.31 and e || = 2.53, and it is 
assumed that the wavelength of the light A = 514 nm. 
Figure ^| shows the order-parameter distribution along 
the boundary for the solution for p = 1. The plot of 
the order-parameter distribution along the boundary is 
an effective way to examine the texture quantitatively. 

When p = 1 , the resulting texture is referred to as the 
boojum texture. It corresponds to a defect with wind- 
ing number +2 [l5| ] lying a distance Rb — l/o^i from 
the center of the domain. As e — ► oo, the virtual de- 
fect retreats to infinity. As e — » 0, corresponding to 
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O 

FIG. 3. The order-parameter distribution along the bound- 
ary shown as a plot of Oo versus tp, where tp is the polar angle 
in the plane-polar coordinates. The parameters are k = 1, 
Ro = 20, and ai = 1.6 

a very strong anisotropic surface energy, or a very large 
domain, the virtual defect approaches the edge of the 
domain. However, the distance of the virtual singular- 
ity from the boundary of a very large domain remains 
nonzero, approaching the limit k/cl\ as e ^ 0. 



For the case of a bubble, instead of Eq. (3.4), we make 
use of 
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Q(k,<p) = 7 [f (e- fe +^) - / (e- fe -^)] , 



(3.8) 



as a solution to Laplace's equation. We find that 

e,(fc, tp) = ) [f a (e" fe+ ^) - f (e- fe -^)] + tt, (3.9) 

(3.10) 



Op 



satisfy Eq. (3.3) in the case of a circular bubble of ra- 
dius Rq centered at the origin. 6^ is show in Fig. |J for 
p = 1, 2, 3, and 4. Also shown in the background 
of each plot is the intensity distribution that would be 
recorded in a BAM image. 

When p = 1, &i in Eq. ( |3.9| ) can be characterized as 
an inverse boojum texture, in that it is obtained by re- 




FIG. 4. The c director distribution and 
flectance for bubble computed for k = 1, Ro = 
where p — 1 in (a), p = 2 in (b), p — 3 in (c) 
(d). 



the BAM re- 
= 20, a p = 1.6, 
, and p = 4 in 



placing k by — k in Qq. This corresponds to a defect 
with winding number —2 located at a distance Rb = oei 
from the center of the bubble. When e = oo, the defect 
is at the origin. As e — > 0, it moves towards the edge 
of the bubble and approaches a distance, n/ai, from the 
boundary. 

Note that in the above discussion the domain and the 
bubble have been assumed to be circular. There is no a 
priori assurance that this shape minimizes the energy of 
the system. 

As the next step, we determine the equilibrium shape 
of the domain or bubble. Rewriting T as p(ip) = e k ^ v \ 
that is, in polar coordinates, we transform Eq. J2.8|) into 



± n b e k + - e)e fc - a"{d - e)e v + [</(tf - e)e v - a"{§ - e)e fc ] k' 

+ W , - e) + s (t - e>] (i - T il 5 ) } + a= o. (3.11) 



where 



d 



tan -1 k' 



for domains 



tp- 



tan 1 kl for bubbles. 



(3.12) 



In Eq. ( |3.11 ) , + applies in the case of a domain while — is 
appropriate to the case of a bubble. Equation (3.11) is a 
nonlinear second-order differential equation, and there is 
no indication that an analytic solution is possible. How- 
ever, for the specific case of domain in which b = and 



a n ^i — except a\, it can be verified that a circular 
boundary centered at the origin is indeed a solution. Fur- 
thermore, such a texture-boundary combination has been 
shown |16| to be a locally stable configuration. Interest- 
ingly, a circular boundar y wit h the inverse boojum tex- 
ture fails to satisfy Eq. ( 3.1l| ) in the case of a bubble. 
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FIG. 5. The geometry of the calculations for domains and 
bubbles in Cartesian coordinates. The gray area is the bulk 
designated by Q. n and t are the outward normal and the 
tangent, respectively. O is the angle made between the c di- 
rector and the x axis and •& is the angle made by the outward 
normal of the boundary and the x axis. 



IV. TEXTURE AND BOUNDARY SHAPE 



In this section, we assume that the virtual boojum 
singularity lies close to the boundary between a domain 
or bubble and the neighboring medium, and we focus 
on the boundary in the immediate neighborhood of the 
singularity. This allows us to treat the two regions as 
semi-infinite. The anisotropic phase is taken to occupy 
(approximately) the half-space for which x is negative. 
The setup of the computation is depicted in Fig. We 
first fix the boundary to lie along the y axis. We then de- 
termine the texture in the anisotropic phase when 6 = 
and all the a„'s except a± are equal to zero. The order- 
parameter field Q(x,y) will be of the form displayed in 
Eq. (3.4), satisfying the boundary condition 
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FIG. 6. The order-parameter distribution and the BAM 
reflectance when the boundary Y is a straight line for re = 1 
and ai = 1.6, where (a) shows Qco for the case of a domain 
while (b) shows Qci for the case of a bubble. 



[kQ x + ai sin(-6)] x=0 = 0. 



(4.1) 



This boundary condition is satisfied by the following ex- 
pression: 



Qco(x,y) = \[h{x + iy) - fo(x- iy)\, 



(4.2) 



where fa(z) is as given in Eq. (3J3) and a = a\/n. This 
texture in fact corresponds to that of a domain in po- 
lar coordinates. Inspection of the boundary condition 
Eq. (4.1) leads us to another solution Qci = — ©Co, 



which corresponds to the texture of a bubble in cylin- 
drical geometry. Figures ||(a) and ||(b) show the distri- 
butions of the order parameter and the computed BAM 
image in this geometry for the domain and the bubble, re- 
spectively. The correspondence between Qco and 6o for 
the case of domains can be observed in Fig. ^(a) in that 
the c directors tend to point towards each other around 
the x axis. Figure ^(b) corresponds to Fig. ||(a) for the 
case of bubbles in that the c directors fan out in the di- 
rection of the outward normal of the boundary near the 
x axis. 

To investigate the equilibrium condition for the bound- 
ary r, we parametrize Y in Cartesian coordinates by 
x = rj{y). The equilibrium condition for T can then be 
written as 



n'Q y ~Q x dd „ 



-|ve| 2 + ( r'(tf-e) _ 

2 T]' 2 at} 



+ ~f7 <T "(fl - e) (e„ + v 'e x ) Vi + ?/ 2 



- [a(0 - 6) + - 6)] 



+ A= 0, 



(4.3) 



and for the boundary condition at y = normal n of T and the x axis. 

, dd , — a 'W lA a\ 

V = -7— ; (l + T) ) — ttt, (4.4) A cusplike singularity occurs when rj (0) ^ 0. As a re- 

sult of the symmetry of the problem, rj(—y) = ?y(y) and 
where •& = — tan -1 ?/ is the angle between the outward rf{—y) = —ffiv)- The possible values of ?/(0) can be 
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obtained by solving Eq. (4.4). That 7/(0) = is a solu- 
tion follows from the fac t tha t cr'(0) = 0. In order that a 
nonzero 77' ( 0) so lves Eq. (4.4), the slope of the right-hand 
side of Eq. (4.4) at the origin must be greater than unity, 
or 



_d_ 

drf 



M 



7(1 + V") 



> 1, 



(4.5) 



which leads us to the cusp condition <r(0) [cr(0) + ct"(0)] < 
0. We will exclude such a condition from our discussion, 
as it requires either |ai| > <7q or 3a2 > <jq. Such con- 
ditions are incompatible with the parameter regime on 
which we focus. 

It can be show n that x = or the y axis is a solu- 
tion to Eq. ( |4.3| ) when 9co is used as the texture (for 
the case of domain) while this is not so for the case of 
the bubble, or when the texture is @ci- We continue 
to investigate perturbatively the response to the texture 
when the boundary deviates slightly fro m x = 0. Let the 
texture, Q(x, y), be of the form of Eq. (3.4) with 



/CaO«/o(*) + £(*)/o(*)- 



(4.6) 



Both t](jj), the deviation of the bou ndary from the 
straight line x = 0, and £(z) in Eq. fl4.6| ) are small quan- 
tities. They possess the following properties: 



r)( x ) = T)(-x) 

e(*) = rco 

The boundary condition can then be expressed as 



(4.7) 
(4.8) 
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d 2 



dx ^ dx 2 



a\ sin(i? — 0) 



(4.9) 



where •& ~ —77'. To first order in 77, £ and their deriva- 
tives, we obtain the following equations: 



2ay£'(iy) = 2a 2 yr](y) 
1 + a 2 y 2 



1 



a 2 y 2 



±0, 



(4.10) 



where -I- corresponds to the case of domain and — cor- 
responds to the case of bubble. The relation between 77 
and £ is readily derived: 



ory ■ 



2ay' 



1 + a 2 y' 2 



±7? 



dy', (4.11) 



Here, we distinguish the primes attached to functions 
which denote derivatives and those attached to variables 
within the integrals which indicate that they are vari- 
ables of integration. Primes will be used in this fashion 
from now on. The function £ obtained from Eq. (4.11) 
is defined along the y axis. We now analytically con- 
tinue £ to the entire x — y plane. Up to this point, we 
have expressed the distortion of the texture in terms of 
a fixed boundary deviation from the exact solution given 
in the be ginn ing of the section. To examine the solu- 
tion Eq. ( 4.11 ), it is necessary to determine the textural 
deformation associated with r)(y) = 770 for domain. We 
find 



£(z) = -arjoz. 
The resulting texture 



Q{x,y) ~ - [fo(x - Tjo + iy) - f (x 
1 



(4.12) 



770 - iy)} (4.13) 



is identical to the domain texture given in Eq. (4.2) ex- 
cept that the position of the virtual defect is translated 
by an amount 770 in the positive x direction. 

Under a small boundary distortion, x — 77 (77) away 



from the y axis, the function f(z) given in Eq. (4.6) and 
its derivative are approximated as 



f(v + ^) = M 1 - ia v) 
f'(v + iy) = - 



1 — iay ' 



1 — iay 1 — iay 
a 2 [77 + £(iy)] 
(1 + iay) 2 



(4.14) 



(4.15) 



where £(i?/) is given in Eq. (4.10). When these ex- 
pressions are substituted into Eq. (4.3) we obtain, for 



© = ©Co, which corresponds to the case of domain, 
2,5a 



V 



1 + a 2 y 2 



[o»? + £'(*»)] 



(4.16) 



while, for 9 = ®a, which corresponds to the case of 
bubbles, 



4a (l - a 2 y 2 ) 
(l + a 2 y 2 ) 2 



2a 2 (3 - 10aV + 3a 2 y 4 ) 8ia 3 y (3 



av; 



(1 + aV) 3 



(1 + aV) 3 



i(iy) 



2a (3 



(1 + a 2 y 2 ) 



(4.17) 



where S = a%/ao. We can see that the right-hand side of 
the equation for the distorting effec t of the texture ap- 
propriate to a domain, i.e., Eq. (4.16), starts a t firs t order 
in 77 while the corresponding equation, Eq. ( 4.17 ), for a 



bubble starts at zeroth order in 77. This provides further 
confirmation that there is no simple inversion symmetry 
between the domain and bubble. 
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V. DOMAINS 



We have established that the boojum texture together 
with a circular boundary is an exact solution for the case 
in which b = and only a\ 7^ 0. This leads us to the con- 
clusion that no nzer o b and/or higher harmonics in the 
expansion Eq. ( 2^2 ) must be present if the boundary is 
to be noncircular. We now attempt to analyze the situa- 
tion in which 01, (J2, and b are all nonzero. We do this by 
perturbing about the boojum solution in terms of small 
parameters b and 7, where 7 = 02/01- 

We note here that the sign of a\ does not affect any 
of the following analysis. As has been mentioned in 
Sec. [0|, if ©o is an equilibrium texture for a\ > 0, then 
= 60 + 7T, representing a reflection of the c directors 
about the x axis, will be the corresponding texture for 
a\ < 0. The equilibrium boundary is circular in both 
cases. Furthermore, contributions of b and ai appear in 
the form of 2 e l2e , which is invariant under reflection of 



the c directors about the x axis. O z represents derivatives 
of O with respect to the variable z which can be x, y, or 
any linear combination of the two. Hence, the effect of 02 
and b is independent of the sign of a%. In the upcoming 
discussions, we will assume a\ > for convenience. The 
inequality 7 > refers to the case in which the c direc- 
tors along the boundary prefer to lie tangent to it while 
7 < applies when the c directors prefer to point along 
the normal to the boundary, ±h. When b > 0, bend 
textures are preferred; splay is preferred when b < 0. 

We hrst find the textural response to 7 and b by mak- 
ing use of Eq. (3.4) with 



f(z) = f (z)+h(z). 



(5.1) 



It can be shown that Eq. (2.4) is satisfied even with 
fi(z) = 0. We continue to investigate the boundary con- 
dition of the texture assuming that the bounding curve 
is a circle of radius Rq. Equation (^) requires a nonzero 
fi(z) satisfying 



K[zffc)-Z- 1 fi(2r)] + % 

nba 
Ro 



1 



z — a 

1 — OLZ 

az 



1 — OLZ 

z — a 

z(l — az* 
z(l - az*) 2 ~ (1 - azf 



hi*) 



a 



1 — az 



1 — az 
z — a 



1 



az 



1 — az 



h(z*) 
2" 



(5.2) 



In contras t to the notation used in Eqs. ( [3.6| ) and (3.7), we have redefined z = e ltp and absorbed Rq into a = R$ai. 
Equation (5.2) can be separated into two equations: 



zf[(z) 



1 



a 



1 — az 



ft(z) = ba 



z(l — az)' 



a 



1 



(5.3) 



and an identical equation with z replaced by its complex conjugate z* . Each of these equations is solvable by standard 
methods. One finds 



h{z) 



7 z — a e(ba — 27) (z — a) 



a(a + e) 1 — az 



2-Fi(l, a/e+ l;a/e + 2; -az), 



(5.4) 



a + e 1 — az 

With fi(z) included, Q(x,y) now satisfies Eq. (|2.4D up to 



where 2F\{v, fi; \x + 1; z) is a hypergeometric function 
first order in 7 and b. 

The full analytical solution of Eq. (3.11) is difficult when 7 and b are both nonzero. However, one can attempt a 
solution as an expansion in the small parameters 7, b, k, and k'. We recall the definition of p{ip) = e k ^ » 1 + k(<p). If 
one ignores terms beyond first order in these quantities, it is possible to solve for the bounding curve, T, analytically. 
The algebraic manipulations are dramatically simplified if we further approximate 6„ » @k/Ro and 6 t ~ 6 v /i?o, 
which is equivalent to neglecting the first-order contributions of k'(tp). The error of the analysis is then of order 
~ 0(Sk'), where 5 = ai/co, which has been defined earlier in Sec. IV. The equation for the boundary T, 



H b R + - e)6 fe - a"(cp - Q)e v 

+ [<r(<p - 6) + a"(cp - 9)] (1 - k")} + A= 0, 



can then be reduced to 



where k" has been separated into various components, as shown below, for convenience, 



(5.5) 
(5.6) 



7 



.... . 5ba 2 (-z + a + e) 

k"(tf) = 7Z - + cc. 



(1 - az) 2 

z 



3 






— a 


2a ~ 


ez) 


(f 


— az 


(r- 


- a 


1 - 


- az 






az 




z 


— a 



fi(z) +C.C. 



(5.7) 
(5.8) 
(5.9) 



ki(tp) consists of terms that depend on 6, /c2(v?) contains terms tha t de pend on 7, and k^{(f) has terms that depend 
on both b and 7 through th e tex tural correction fx given in Eq. ( |5.4| ). The functions k'(ip) and k(ip) can then be 
obtained by integrating Eq. (5.6) with respect to ip. They are 



k[(ip) 
k' 2 {ip) 



2 r 



Sba 

i 

5a"f 
i 



— (e + a) ln(l — az) + 



2a + 2 



1-a 4 



1 — a — ae 
a[az — 1) 

ln(l — az) 



cc. 



J_ 3a 1 - a 2 

2^ ~ T + e 



1-a 2 (1-a 2 )' 
az — 1 2a(az — l) 2 



cc. 



5a 



2l 



where 



i F 1 (l,a/e+l;a/e + 2;-z'), , 
: az . 



One more integration yields 



Sba 2 


[( 


1 

e 






a 


Sj( 




'2 

- + 


a \ 




a 


-a(l 


-a 2 ) 



+ a I ln(l — az) — (e + a)Li2(az) 



cc 



3_a^ 
~2~ 



1 - a 4 



Li2(az) + 



(l-a») S 



ln(l — az) — 



1-a 2 



2a(az- 1) 



cc, 



2l 



kziv) = [ba — h(az) + cc, 



where 



Hz) 



g(z')/z'dz' 



and Li2(z) is the polylogarithmic function defined as 

Li„(x)=^^ 



OO U 



k=l 



Numerical integrations can be utilized for the evaluation of k(ip). However, as we can see from the following, 

z z 'a/e „a/e+l 

dz = 2-Fi(lj a/e+1; a/e + 2; — az), 



1 — az' 



a + e 



(5.10) 

(5.11) 
(5.12) 

(5.13) 

(5.14) 

(5.15) 
(5.16) 

(5.17) 
(5.18) 

(5.19) 



the integrand oscillates strongly as a/e — ► 00 as a result of the factor z' Q / e and numerical integrations become 
inefficient. Further observation reveals that g(z) and h(z) can be evaluated analytically if a/e = n is an integer. 
Equation (5.19) simplifies as follows: 
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1 — az 1 



-dz' = -a""- 1 



E 



{az) 1 



ln(l — az) 



(5.20) 



When the above simplification is substituted into Eq. ( 5.13| ), the integration can be performed and yields the analytic 
form of g(z): 



/(*) = E 



n — k + 1 n+1 



1 ln(l - z) 
kz k n + 1 



yU+1 



(5.21) 



h(z) can be evaluated analytically in the same manner and we get 



Kz) = 



(n — k + l)k 2 z k 



E 



i 



i 



(n+l) 2 k (n+l)k 2 



ln(l-z) ln(l-z) Li 2 (z) 



(n + l) 2 z n+1 (n + 1) 2 n+1 



(5.22) 
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The boundary of the domain fc(y) = / y fcj (y) is 

smooth and has continuous derivatives. Typical fc(y), 
k'(cp), and k"(ip)'s are shown in Fig. 0(a). The corre- 
sponding boundary T parametrized by p(<p) = 1 + k(ip) 
is depicted in Fig. |?](b). We have thus arrived at an 
approximate expression for T as a function of the line- 
tension anisotropy coefficient 7 and the elastic anisotropy 
coefficient b. This expression is useful when we are inter- 
ested in the response of T for small values of the these 
anisotropic parameters. 

We first examine the boundary response to 7 while 
keeping 6 = 0. We find indentations and protruding fea- 
tures on the domain boundary for 7 < and 7 > 0, 
respectively. The progressive change of the boundary re- 
sponse when ci2 changes from —0.5 to 0.5 is illustrated in 
Fig. ||. The results are in qualitative agreement with 
those presented in Ref. ]lO|j . We have also examined 
the dependence of the boundary at fixed 7 on domain 
size, i?o- Figure ^ show boundaries for domains of sizes 
i?o = 0.2 to 10. When R < 1, the domains appear 
slightly flattened and elongated if 7 > 0. This is in ac- 
cord with the intuitive notion that the second-harmonic 
term in the line tension becomes important as the vari- 
ation in the texture vanishes, i.e., in the limit that the 
order parameter is uniform. As the domain gets larger, 
the boojum singularity moves closer to the edge of the 
domain and the boundary correction moves towards the 
axis connecting the center of the domain and the boo- 
jum. Cusplike features start to appear when the domain 
is larger than a "threshold" size, i?o = 1 for the domains 
in Fig. ^ We note that we have used large values of 
7 < 0.3 to illustrate the nontrivial boundary that we 
have obtained for the domains. It has been numerically 
verified |l3| that the qualitative behavior of the bound- 
ary response is indeed preserved up to much larger values 
of 7. 

It has been shown in Sec. |y| that the boundaries of 



the domains are strictly smooth and continuous in the 
parameter regime upon which we focus. We are, how- 
ever, able to find domains with cusplike features in the 
context of the perturbative analysis described in this pa- 
per. Such domains can be characterized by an excluded 
angle *S>o defined in Fig. |To|(a). Domains with boundaries 
that resemble those with cusplike features are observed 
experimentally. The domain-size dependence of "to is 
shown in Fig. |l0|(b) g. There is no rigorous mathe- 
matical definition of ^0 for a continuous boundary. It 
is, nevertheless, possible to devise a systematic way of 
identifying such an excluded angle for a smooth bound- 
ary. One first evaluates W = —2 tan -1 dx/dy along the 
boundary. The value of ^ at the straightest part of the 
boundary, which is indicated by d 2 x/dy 2 — > 0, is a likely 
candidate for Wo- Figures |ll|(a) and |ll|(b) show the plot 
of '5 and d 2 x/dy 2 versus ip, respectively. The values of 
"J in the plateau region in Fig. |ll](a) represent the range 
in which the measured excluded angles are likely to fall. 
These values of 'J are found in the region near the $ axis 
of the plot of \t versus d 2 x/dy 2 shown in Fig. |TT](c) . A 
plot of '5 versus / = Iq exp[— (d 2 x/dy 2 ) 2 ], as shown in 
Fig. |ll|(d), highlights the range of the values of ^ that 
is most likely to contain the measured excluded angle. 
Figure |ll|(e) displays / as the intensity (inverted, in that 
the brightest corresponds to the smallest value of I). Fig- 
ure |ll](f) shows the value of ^ at which / = Z m ax by 
the dark line and the region in which / > /max/2, or the 
full- width- at-half-maximum (FWHM), by the gray band. 
Figures |lT](e) and [ll](f) are useful for describing the se- 
lection process by which one is led to the most likely 
values of the excluded angle. Figures |l|(a) and |l2|(b) 
illustrate such plots. The experimental result is super- 
posed in Fig. |l2|(b). The parameters are adjusted in order 
to obtain a by-eye fit. The values of the parameters are 
nja\ = 4 um, 6 — 0.4, and 7 = 0.5 

The perturbative analysis generates results that arc in 
good agreement with experimental observations |J for 
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FIG. 7. (a) Plot of k"{tp), k'(ip), and k(tp) for « = 1, 
i?o — 5, (To = 4, ai = 1.6, and a2=0.5. (b) The corresponding 
domain shape T parametrized as p(ip) « 1 + 

large domains. It has also captured qualitatively the es- 
sential features, namely the onset and the maximum of 
the domain size dependence of the excluded angle. As 
displayed in Fig. [13(b), the maximum and the onset of 
are quantitatively different in the perturbative anal- 
ysis and in the experimental data in the intermediate Rq 
regime. In particular, the experimental maximum of \I/o 
cannot be obtained from the analysis even though 7 = 0.5 
has been used. We shall defer the discussion on this to 
the end of this section after elaborating on the effect of 
the elastic anisotropy b. We also note that 7 = 0.5 is very 
large as a perturbative parameter. Although there is no 
a priori guarantee of the accuracy of the results, it is evi- 
dent from our numerical studies jl3| that the qualitative 
behavior of the boundary as a function of the domain 
size is preserved in the perturbative analysis up to at 
least 7 = 0.5. 

We now proceed to examine the effect of b on the 
boundary T. The coefficient of the anisotropic line ten- 
sion 7 is kept at zero and the boundary response is pro- 
portional to b. Figure [l3] shows the plot of k(ip)/b. In 




FIG. 8. Domain shapes computed for k = 1, Rq = 5, 
ao = 4, ai = 1.6. and (a) a 2 = —0.5, (b) ai = —0.3, (c) 
a 2 = -0.1, (d) a 2 = 0.1, (e) a 2 = 0.3, and (f) a 2 = 0.5. 




FIG. 9. Domain shapes computed for k = 1, ao = 4, 
ai = 1.6, a 2 = 0.6, and (a) R = 0.2,0.25,0.33,0.5,1 and 
(b) R = 2,2.5,3.3,5,10. 

contrast to the results obtained by Galatola and Fournier 
||10|| , the boundary acquires a denting correction when 
b < 0, indicated by a maximum in k(ip)/b at tp = 0. 
This perturbative result is confirmed for small values of 
b (=0.1) by our numerical studies p8[ . At such small 
values of b, the boundary is practically circular. We shall 
restrict our discussion to the effect on T of small values 
of b, in that higher-order corrections of b, which are not 
taken into account in this first-order perturbative anal- 
ysis, change the qualitative behavior of the boundary 
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R ((xm) 




l/R (l/|xm) 

FIG. 10. (a) The definition of the excluded angle * - (b) 
Experimental measurements of the domain-size dependence 
of observed in L2 domains surrounded by LE phase taken 
from Ref. [5]. 
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FIG. 12. (a) Density plot of J as a function of '3/ and Ro. 
(b) Plot of /max and the region in which / > /max/2 as 
a function of ^ and Rp. Superimposed are the experimen- 
tal data shown in Fig. Iicj(b) with parameters /t/ai =4 fim, 
5 = 0.4, and 7 = 0.5. 
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FIG. 11. (a) Plot of * = -2 tan" 1 dx/dy versus ip for 
k = 1, Ro = 5, ao — 4, ai = 1.6, and 0,2 = 0.6. (b) Plot 
of d 2 x/dy 2 versus tp for the same parameter, (c) Plot of ^ 
versus d 2 x/dy 2 . (d) Plot of ^ versus / = exp[— (d 2 x/dy 2 ) 2 ]. 
(e) Density plot of / as \l/ for a single _Ro- (f) The dark line 
marks the maximum /max of / while the gray region shows 
the range of in which / > /max/2. 

response, as reflected by our numerical studies |13l. 



0.03 




-0.01 

FIG. 13. Plot of k(ip) /b for k = 1, Ro = 5, a = 4, 01 = 1.6, 
and <i2 = 0. The maximum of k(tp) at ip — implies a pro- 
truding correction when b > 0. 

As for the role of b in the interpretation of the ex- 
perimental observations, we conclude in our numerical 
studies [ll| that a nonzero value of b cannot be solely 
responsible for the protruding features, and hence the 
excluded-angle measurements. Large values of 7 ~ 0.5 
are required to produce excluded angles whose maximum 
approaches the largest value of experimentally measured 
excluded angles. At sufficiently large values of 7, we have 
found that the behavior of the excluded angle is quali- 
tatively unchanged when the anisotropy parameter b is 
varied from —0.5 to 0.5. Although the validity of the 
perturbative analysis at b > 0.1 is questionable, the rela- 
tive magnitude of the correction to the boundary of b is 
much smaller than that of 7. This is further verified by 
numerical studies ||i"8|| . 

We have thus demonstrated, within our first-order per- 
turbative analysis, that the line-tension anisotropy 7 can 
give rise to the indentations and protruding features of 
the domain boundaries that have been experimentally 
observed |Q. Our results on the boundary response to 
7 are in qualitative agreement with prior results Jl0| , ^3| . 
Although our investigation of the b dependence of the 
boundary does not provide us with dependable results 
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for large values of b, it supports the assertion that the 
textural correction is an important contribution to the 
boundary response. Further analysis of the available ex- 
perimental data on the excluded angles leads to the con- 
clusion that b has little effect on the boundary of the 
domains with protrusions. We shall confine our conclu- 
sions to small values of \b\ < 0.1, although large values 
of b ~ 0.8 do lead to interesting domain shapes. Discus- 
sions of the domain boundaries at large values of b will 
be presented in a forthcoming article |l8[ , and some of 
the results have been briefly presented in Ref. |l3| . 

We obtain good agreement between the results of the 
perturbative analysis and the experimental observations 
on the dependence on the domain radius Ro m the 
large- Rq regime. In the intermediate- Rq regime, the dis- 
crepancy is not resolved, even in our numerical stud- 
ies ||l8| . The mismatch could possibly be attributed to 
the fact that our simple elastic theory does not describe 
accurately the actual complex monolayer. In the per- 
turbative analysis performed in this work, the parame- 
ters are restricted to a region in which a(ip) + <j"(ip) is 
always greater than zero. Furthermore, it is generally 
known that the dipolar interactions between the surfac- 
tant molecules in the monolayer are important. The cur- 
rent model does not take into account such interactions. 
It has been discovered in a recent experimental study fl9| l 
that the tilt is not always uniform, especially in the re- 
gion around a point defect. The contribution of variation 
in tilt may not be significant in terms of accounting for 
the discrepancy in the intermediate-sized domain regime. 
It does become important in the large- Rq regime when 
the virtual singularity approaches the domain boundary 
and the texture acquires a rapid variation in the neigh- 
borhood of the virtual singularity. 



VI. BUBBLES 




(bl 



FIG. 14. Bubble shapes computed for k = 0.16, ao = 1, 
«i = 0.16, and (a) R = 1, 2, 3, 4, 5, 6, 8 and (b) R = 8, 20, 40. 
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FIG. 15. (a) Plot of / as a function of * and Ro. (b) Plot 
of /max and the region in which I > /max/2 as a function of 
$ and Ro. Superimposed are the experimental observations 
of gaseous bubbles in L2 phase. The experimental data have 
appeared in Ref. [5]. The parameters for the by-eye fit are 
n/a\ — 0.4 /jm and 5 — 0.16. 



26e 
i 



1 



+ ln(l - az) 



1 — az 

k((f) = —2Se [ln(l — az) + Li2(az)] + c.c. 



(6.2) 
(6.3) 



It has been shown in Sees. Ill and IV that there is 



no straightforward inversion symmetry between the do- 
mains and the bubbles. In contrast to the case of the 
domain, it is not necessary to introduce anisotropic pa- 
rameters other than a\. The inverse boojum texture ®i 
given in Eq. ( |3.9| ) and Eq. ( |3. 10 ) satisfies the equilibrium 
condition for a circular bubble for the case in which the c 
directors favor pointing into the bulk or a\ > 0. Substi- 
tuting Q; in to the equilibrium condition for the boundary 
r Eq. ( |3.11 ) , one finds that a circular boundary does not 
satisfy the equation. By perturbing about the circular 
boundary in terms of a small parameter 5 = ai/co, we 
arrive at an equation similar to Eq. (5.6), 



k"(<p) = ~2a 2 Se 



(1 — az)' 



+ c.c. 



(6.1) 



Again — see Eqs ( |3.9| ) and ( |3.10 ) — we have redefined z = 
e lv and a = ai/Ro. Following the same procedure as for 
the case of a domain, we find for fc'(^) and k(ip) 



It should be kept in mind that the above discussion refers 
to the case a\ > 0. The results for a\ < are not ob- 



tained by a simple sign reversal of a\ in Eq. (S.3). Appro- 
priate changes in the texture and definition of e, which 
is given as k/(\cli\Rq), must be taken into account. The 
details of the calculations are presented in Appendix |^. 
The final bounding curve for the bubble depends only on 
the magnitude |ai|. We have derived expressions for T 
for the cases where there is only an a\ contribution in 
the line tension. 

One can utilize the results to investigate the bubble- 
size dependence of the boundary shapes. Figure [l^ shows 
the shapes for the bubbles of sizes Rq = 0.2 to 10. Very 
small bubbles i? <C 1 are nearly circular, as are very 
small domains. Cusplike features start to appear when 
the bubble is larger than a "threshold" size, Rq — 1 for 
the bubbles shown in Fig. Similar analysis of the ex- 
cluded angle to that for the domain presented in Sec. |V] 
can be carried out. Figure shows plots of ^ versus Rq 
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corresponding to those in Fig. [lj. Experimental mea- 
surements H of the cusp angle for the bubble are super- 
posed in Fig. [l5|(b). A by-eye fit can be obtained with 
parameters n/ai = 0.4 /mi and S — 0.16. We find fairly 
reasonable agreement between the theory and experimen- 
tal observations. We remark that the apparent mismatch 
between the theory and the experimental data point at 
i?o = 1 A 4111 ; which has been shown to match the explicit 
measurement on the calculated bubble boundary in Ref. 
|^| , may be the result of the inadequacy in qualifying the 
excluded angle for bubbles with R < 1 /zm using the 
FWHM of I shown in Fig. ||(f). 

As compared to the parameters obtained for the do- 
mains in Sec. |y|, which are k/cl\ — A /im, S = 0.4, and 
7 = 0.5, the value of n/ai for the case of bubbles is an 
order of magnitude smaller than that for the case of do- 
mains. Noting the fact that the data for the bubbles 
are obtained at the L 2 /G coexistence region and those 
of the domains are measured when the L 2 domains are 
surrounded by the LE phase !|, the comparison of n/ai 
between the domains and the bubbles is indeed in accord 
with the intuitive sense that n of the L 2 domains should 
not vary significantly while a\ at the L 2 /G interface is 
much larger than that at the L 2 /LE boundary. The <5's 
are of the same order of magnitude and there is no cor- 
responding 7 in the case of bubbles. The result of the 
perturbative analysis is consistent between the domains 
and the bubbles. 



fluctuations and its implications to the computation that 
has been carried out. The effect of fluctuations can be 
assessed by utilizing a mapping between the statistical 
mechanics of the order-parameter fluctuations in this sys- 
tem and the behavior of a two-dimensional Coulom b ga s 
pof . Consider the Hamiltonian of the form of Eq. (|2.l|), 
with 6 = and a n — for all n ^= p. For a system with 
circular boundary of radius Rq, one has 



H'[&\ = H[Q] - 2na Ro 



(7.1) 



the prime here distinguishes the free energy from the one 
defined in Eq. ( |2.1| ). The prime is dropped from now on 
for convenience. The partition function can be written 
as 



Z(a p ) = Z(0) ^exp 



(3a p <£ cos p(ip — Q)ds 



(7.2) 



where (O)o denotes the thermal average with respect to 
the Hamiltonian without the boundary term given below, 



(O)o 



/I?eOexp(-^ / dA\VO\ 2 ) 

/peexp(-^/dA|ve| 2 ) ' 



(7.3) 



We denote 0&(</?) = @(#o costp, i?o simp) as the val- 
ues of on the boundary of the circular domain. The 
following correlation function can be evaluated |Ll| , 



VII. THERMAL FLUCTUATIONS 



-A In 2 



tp-ip 



(7.4) 



The analysis presented in the earlier sections is in the 
mean-field approximation; thermal fluctuations are ig- 
nored. In this section, we examine the effect of thermal 



where A = 1/2it[3k. 

To evaluate the full partition function Eq. (7.2), we 
Taylor expand the exponent as 



Z{a p ) 
Z(0) 



71=0 ' \i = l 



(7.5) 



where we have added an index i to distinguish the various cosine terms, denoted 0^ = Q(Rq cos ipi, Rq sin ipi) for 
convenience, and introduced a microscopic length scale a in the denominator in ds to make the total integration 



dimensionless. We use cos 1 

Zjflp) 
Z(0) = 



(« 

E 

n=0 



,i<t> 



-i<f>\ 



-/3% 



/2 and then expand the products of the cosine terms, 



(7.6) 



where we have defined charge qt = ±p for each of the cosine terms in Eq. (7.5) and denoted a sum over all charge 
configurations. The thermal average can be evaluated exactly using Wick's theorem; we obtain the following equality: 



(7.7) 



The average (0 2 )o encountered in Eq. (7.7) can be evaluated using the inverse of the microscopic length scale 1/a as 
the ultraviolet cutoff. We have the following: 
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(e 2 M)o = um (e(^)e(^))o. 

tp-npi 

= Aln — . (7.8) 
a 

The contributions of the charge configurations in which 7^ to the partition function are suppressed as a result 

In — I . Therefore, it is necessary to sum only over the configurations in which there 

2 a '. . 

is no net charge. Finally, the full partition function can be expressed as follows: 



analogous to the expression of the partition function for a 
system of 2 D neu tral coulomb gas. The charge of the par- 
ticles is py/ A/2 and the particles are distributed along 
the circumference of the domain. Following the Coulomb 
gas treatment for the 2D phase transitions |H| , we derive 
flow equations for the running coupling constants, 



da p (a) 

da 
dA(a) 

da 



1 - 



p 2 A(a) 



o P (a), 



-2tt 2 A 2 (a)a 2 p {a). 

For a p (a) < 1, we have the following: 

1 \ / \ i-p 2 a/2 
a p {a) - " x 



(7.10) 
(7.11) 

, , v , • ( 7 - 12 ) 

The scaling relation implies the relevancy, and the irrel- 
evancy transition temperature of a p (a) at T p is given by 

k B T p = — r . 7.13 

As compared to the Kosterlitz-Thouless transition tem- 
perature ksTxT = 7tk/2, T\ and T2 are above Tkt- Our 



result is analogous to the scaling index of the symmetry- 
breaking perturbation in the 2D planar model obtained 
by the spin- wave approximation [ p2[ . 

In the low-temperature phase, we consider fluctuations 
up to a cutoff that is proportional to the sample size, or 
a ~ Rq, and consider the renormalized coupling constant 
a p = a p (Ro)Ra- We find the scaling relation, 



R 



-p 2 A/2 



(7.14) 



Based on a theory of fixed K, the renormalized anisotropic 
line tension decreases as a power law with Rq, the radius 
of the boundary, with exponent — p 2 A/2. Using this re- 
lation, we investigated the effect of thermal fluctuation 
on the domain boundary. The result is depicted in plot 
of ^ at the maximum of /, VPo, as a function of Rq in 
Fig. |l6]. We notice the rounding off at the maximum of 
^0 and the decrease in the magnitude of ^0 with increas- 
ing temperature. The above comparisons are made for 
domains with S = 0.4 and 7 = 0.38 when Rq — 0.2. 

We can also look at the effect of fluctuations on b. We 
let a n — for all n and the Hamiltonian becomes 



H[Q] = ^ J { I V9| 2 + b [(-6 2 + 9 2 ) cos 26 - 26^ sin 26] } dA. 



The partition function 



Z(b) = Z(0) <^exp 



f3bn 



J (6 2 e i2e + 6 2 e- 42e ) dA 



(7.15) 



(7.16) 



which leads to the flow equations similar to those for a^, 
db(a) 



da 
dA{a) 
da 



-2A(a)b(a), 
-27r 2 A 2 (a)6 2 (a), 



(7.17) 
(7.18) 



from which we deduce that b is always irrelevant at finite 
temperature. 



Following the same argument as for the case of a p , 
we consider fluctuations that cut off at a ~ Rq in the 
ordered phase, and take the renormalized coupling con- 
stant 6 = b(Ro). We find that the b scales as Rq 2A in 
a theory with fixed n. When thermal fluctuations are 
important, the magnitude \b\ is at its maximun (< 1) for 
the smallest Rq and decreases as Rq increases. Thermal 
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FIG. 16. Plots of as a function of Rq for « = 1, ao = 4 
at T = 0, O.ITa-t, O.STkt, 0.5T kt , and 0.7T KT . At i? = 0.2, 
all domains have the same renormalized coefficients in the ex- 
pansion of the anisotropic line tension, namely ax = 1.6 and 
a 2 = 0.6. 

effects reduce the already insignificant boundary correc- 
tion due to b. This reinforces our earlier conclusion that 
b has no significant influence on the domain boundary. 

In summary, thermal fluctuations act to renormal- 
ize the anisotropic parameters. The influence of ther- 
mal fluctuations on the boundary shape can be stud- 
ied using the mean-field approximation with renormal- 
ized anisotropic parameters. The boundary correction 
due to elastic anisotropy, which is negligible at T = 0, 
further decreases as a result of thermal fluctuations. 
The deviation of the boundary from a circle that re- 
sults from line-tension anisotropy is displayed in terms 
of the domain-size dependence of the excluded angle in 
Fig. |h]. The maximum of the excluded angle, which is 
higher in the experimental observations than predicted 
at T — (see Fig. |l2]), will be reduced when thermal 
fluctuations are introduced. This indicates that the line- 
tension anisotropy of the monolayer system under study 
may be very strong. 



VIII. CONCLUSIONS 

We have described in this paper a systematic investiga- 
tion of a system of a 2D ordered mediu m wi th a bound- 
ary. Beginning with the free energy Eq. (2.1), which de- 
scribes a bounded monolayer, we have derived the Euler- 
Lagrange equations for the texture and the boundary. 
From the boundary conditions we have shown that the 
boundary under consideration must be smooth. A con- 
tinuation of the analysis in the spirit of Ref. Q for the 
cases of both the domains and the bubbles reveals that 
bubbles do not remain circular when the only term in 
the anisotropic line tension, as given by Eq. ( |2.3| ), is 
ai cos 4>, while circular domains are not affected. There 
is, thus, no simple inversion symmetry between domains 
and bubbles, as one would intuitively expect. Pertur- 
bative calculations have been carried out to investigate 
the influence of the term, the second-harmonic con- 



tribution in the line tension, and the elastic anisotropy, 
parametrized by b, for the domains. Assuming the do- 
mains are nearly circular and the anisotropics are weak, 
the perturbed domain shapes are then computed analyt- 
ically to first order in small parameters. We have exam- 
ined the boundary response to the anisotropic parame- 
ters 6 and a.2- Our results for the boundary response to 
a.2 are in qualitative agreement with those reported in 
Refs. pcj , p^ |. We have also obtained the dependence of 
the boundary shape on domain size that is in qualitative 
agreement with experimental findings when 6 = 
As for the boundary response to 6, our perturbative re- 
sults contrast with the conclusions arrived at in Ref. [|l0| . 
Textural correction plays an important role. These con- 
clusions are confirmed in our numerical studies |]l3| , |l8| . 
The b contribution to the boundary is much weaker than 
that of the line-tension anisotropy. The quality of the 
fit to the currently available experimental observations 
is not sensitive to the value of b in the range from —0.5 
to 0.5. Considering only the line-tension anisotropy 02, 
the result of the perturbative analysis has qualitatively 
captured the essential features of the experimental obser- 
vations. The quantitative mismatch can be attributed in 
part to the inapplicability of a perturbative approach to a 
parameter region in which the anisotropic parameters are 
large. The detailed difference between the simple model 
we adopt and the actual complex underlying structure of 
the L2 domains may also contribute to such discrepan- 
cies. Long-range dipolar repulsion has been ignored. The 
tilt degree of freedom , which may not be significant 
in the small domain regime but can be important in the 
large- Rq region, is not included. 

In the case of the bubbles, we evaluate the boundary 
response due to the a\ contribution in the line tension. 
We are able to produce a dependence of the shape of the 
bounding curve on bubble size that compares favorably 
with experimental observations. This result has been re- 
ported earlier in Ref. Q . The parameters of the domains 
and those of the bubbles, both obtained with by-eye fit of 
the result of the analysis to the experimental data, are in 
reasonable agreement, taking into consideration the fact 
that the data for the domains are taken in the L2/LE 
coexistence region and those for the bubbles is obtained 
for gaseous bubbles surrounded by the L2 phase. Fi- 
nally, we have presented an analysis of the effect of ther- 
mal fluctuations that leads to a domain-size dependence 
of the line-tension and elastic anisotropics. The analy- 
sis also suggests that the line-tension anisotropy may be 
very strong. 
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APPENDIX A: GEOMETRY AND COORDINATE 
SYSTEMS 

We enumerate in this section the forms taken by vari- 
ous geometrical quantities of a 2D space curve in different 
coordinate systems and the relationships between those 
forms. A curve T surrounding £1 can be represented by 
a one-parameter trajectory of the position vector r(i), 
where t is parameter. Its unit tangent vector is given 
by t — dr/ds and the unit outward normal is given by 
n = ri/\ri\, where n = ±d 2 r/ds 2 and ds = \df/dt\dt. 
We let i) be the angle between the normal vector of the 
curve and a reference axis. Then the radius of curvature 
is \ds/d&\. 

Consider the problem of a nearly circular domain. The 
coordinate system of choice is obviously plane-polar. It 
is convenient to use tp as the independent variable. We 
then write T as f(<p) = e fc ^ p, where & k ^ v > is the radial 
distance from the origin and tp is the polar angle at r. 
A length element given by ds = e fc Vl + k' 2 dtp is in the 
positive direction of f2. The unit vector n = —d 2 r/ds 2 
points away from the origin, or outwards from T, 



k'p + tp 

VI + k 12 

p — k' 'tp 

VI + k' 2 ' 



(Al) 
(A2) 



We have cost) = n 
follows that 



d0_ 

ds 



x, which gives d = tp — tan 1 k' . It 



k" 



1 + fc'V e fc VTT^ 



(A3) 



Similar relations can be derived for the case of a nearly 
circular bubble. Here, the length element in the positive 
direction of U is ds = — e k Vl + k' 2 dtp and the outward 
normal n = d 2 r/ds 2 . We obtain d = tt + tp — tan -1 k' . 
The geometrical quantities are evaluated as follows: 



dO 
ds 



=^ = - 1- 



k'p + tp 

VI + k' 2 ' 

p — k' tp 

Vl + k' 2 ' 
k" 



1 



l + k' 2 J e fe Vl + fc' 2 ' 



(A4) 
(A5) 
(A6) 



Cartesian coordinates are useful when we are interested 
in a small region of a large circular domain or bubble, on 



the scale of which the boundary is nearly a straight line. 
In Cartesian coordinates, we have for the position vector 
r(i) = x(t)x + y(t)y, where x and y are the unit basis vec- 
tors. We can always pick y as the independent variable 
and x = rj{y) as the dependent variable for the curve T. 
The symbol rj(y) is chosen deliberately to avoid confusion 
with the independent variable x for the texture O. We 
take f2 to reside in the region x < 0, assume that T nearly 
coincides with the y axis, and take the virtual defect to 
be in x > near the y axis as depicted in Fig. ||. We have 
ds = yl + rj l2 dy traversing along the positive direction 
of r and n = —d 2 r/ds 2 . We immediately obtain 



t = 



T) x + y 
x - rj'y 



V 1 + n 
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(A7) 
(A8) 



The angle between the normal and the x axis is $ 
— tan -1 77', and 



M 
ds 



V 1 + i' 



(A9) 



APPENDIX B: THE EXTREMUM EQUATIONS 



We begin with the free energy Eq. (2.1). In Ca rtesi an 
coordinates, the elastic energy density Hb in Eq. (2.2) is 
given as 



H b = ^{\ve\ 2 + b[(-e 2 x 



e;) cos 2e- 



-e x e y sm2e}} . (bi) 

Taking the variation of H[Q] with respect to 0, we find 
~dH b d dU b d &H b ~ 



5H = 



dS dx dO x dy dO y 
dHb dx dHb dy 



(59 dA± 



d<d y ds 



dO y ds 



SO ds~ 



</(??- 6)£6 ds. 



(B2) 



The boun dar y integral <fL is taken counterclockwise. The 



+ in Eq. (B2) is appropriate to the case of domains while 
— is appropriate to bubbles. The equilibrium condition 
SH/80 = results in the Euler-Lagrange equation 



d dlib d dm 



dHb 



86 dx de x dy d@ y 







(B3) 



for (x,y) € f2, which can be reduced to Eq. (2.4). The 
boundary conditions can be expressed in Cartesian coor- 
dinates system as follows: 



1G 



f _ dy ~ dx 



e — + e ^ 

<is ds 



cos 26 



e — -e ^ 

as as 



sin 26 



± </(i? - e)= o, 



(B.4) 



To display the derivation of the equilibrium equation for the bounding curve T, it is more convenient to utilize 
plane-polar coordinates. In the case of a domain, we rewrite the free energy in plane-polar coordinate as follows 



/7T rk 
/ H 6 e 2fel dfci + a(ip - tan" 1 k' - 0)e k \/l + 
-7T JO 



(B.5) 



We then take a variation of the free energy with respect 
to k((f>). The equilibrium condition results in the Euler- 
Lagrange equation 

OH d &H „ 

(B.6) 



5fc aV dk' 



0. 



where 



Jo 



a(<p - tan" 1 k' - 9)cVl + k' 2 . (B.7) 



aw 

9fc 



H b e 2k + {a-a'Q k )e k ^/l + k' 2 . (B.8) 



We have also the partial derivative of 7i with respect to 



dU _ (-a' + ak')e k 

dk' ~ VTTfc 72 " ' 



(B.9) 



To continue, we now look at the partial derivative of H Taking the derivative of Eq. flB.9|) with respect to the 



with respect to k. We have 



independent variable tp results in 



d on 

dip dk' 



a" (6 fc fc' + e v ) - a'k' (Q k k' + %) ak' 2 + a" {a + a")k" 



VI + k' 2 



vTTfc 72 " yTTfc 72 " 3 



(B.10) 



Equation (3.11) for the case of a do main follo ws imm edi- 
ately from the substitution of Eqs. (B.8) and ( B.10 ) into 



the Euler-Lagrange equation, Eq. (B.6). We have as- 
sumed here that the curve joins smoothly at (p — ~tt,tt 
and there is no boundary contribution from these two 
end points. However, we are particularly interested in 
finding out if a cusp, in the form of a discontinuity in 
the slope of the bounding curve, exists. In our system, 
which is symmetric about the x axis, the singularity is 
expected to occur on the x axis. We thus allow for the 
possibility that T has a discontinuity in slope at ip = 
and determine the condition for the discontinuity. The 
assumption of a discontinuity in T gives rise to an extra 
boundary condition at ip = 0, 



OH 
dk 1 



0. 



(B.ll) 



Using Eq. (B.9) and the fact that &(e k ) = 0, we thus get 
er't-tan™ 1 *;') 



cr(- tan -1 k') 



(B.12) 



For the line tension given in the fo rm of Eq. (2.3), k' = 
is always a solution to Eq. (B.12). In order for k' to be 
nonzero at <p = 0, it is necessary that the slope of the 
right-hand side as a function of k' at k! = be greater 
that unity, i.e., 



d a'(-tan- 1 fc') 



dk' cr(- tan" 1 k'^ 



> 1. 



(B.13) 



k'=0 



This implies cr(0) [ct(0) + a"(0)} < 0. We have excluded 
such a condition, in that it requires either |ai| > ctq or 
3a2 > (To, both of which are well beyond the parameter 
regime that we are focusing on. The boundary T that we 
are solving for will not have a singularity. 
In the case of a bubble, we have 



H[<p; k, k'] 



Hte^dh 



a(n + ip- tan" 1 k' - 9)eVl + k' 2 . (B.14) 



This is similar to Eq. ( B.7 ) apart from the limit of in- 
tegration in the first term on the right-hand side of the 
equation. 

In the case of the Cartesian coordinate system, we 
write the free energy as 



H[ V ] 



H b dx + ct(i? - 6)v/l + V 2 



dy (B.15) 



and the equilibrium condition can be obtained using the 
following Euler-Lagrange equation: 
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dU _ d_&H _ 
drj dy drj' 



(B.16) H((p\ri,rf)= p H b dx + a(-d-e) v / l + T 1 ' 2 . (B.17) 

J — oo 



where 



APPENDIX C: SAMPLE CALCULATION 



We present here an analysis of the equilibrium equations for the special case in which b = a nd only a single a p > 0. 
In this case, the bulk equation is automatically satisfied if we write 6 in the form of Eq. (^4). We first consider the 
case of a domain, for which the boundary condition is given by Eq. ( |3.2| ). We substitute the boojum texture Oo and 
find 



KCXC 



pk 



iR 



-ip^p 



K 

iR 



1 - aeP( fe+i v) 
1 



1 - aeP( k+ll f) 
(a 2 e 2pk + 2eae pa - 1) 



1 — aeP( fe_i ^) 
1 

1 - ae P(k-i<p) 
1 



a p i 
~2~ 



2aeP k 
1 



,J>k 



-iptp 



„pk 



1 - aeP { - k+i ^ 
1 



a 2 e 2 P k 



1 - aeP( k+i ^ 1 - aeP^-W 



1 - aeP( k+il P) 
0. 



1 - aeP( fe -^) 
1 - a 2 e 2 P k 
1 - aeP^-^ 



0, 

= o, 



(C.l) 



The above boundary condition is to be satisfied for all (e , tp) € T and hence the coefficient in Eq. (C.l) has to be 
0, which gives ae pk = — e ± yl + e 2 , w here e = k/ (pa p Ro). Together with the requirement that © does not have a 
singularity in 17, we arrive at Eq. (3.7), where Rq = e k is the radius of the circular boundary. As for the case where 
a p < 0, we substit ute Q_ = ©o + ir/p into the boundary condition Eq. (3^). Requiring that the texture is continuous 
in f2, we find Eq. (3.7) with e = K,/(p\a p \Ro). The results for bubbles can be obtained in a similar manner. Although 
we have demonstrated a solution of O for p any integer, a circular domain shape does not satisfy the equilibrium 
condition for T in general. 

We now examine the domain shape for the case in which p = 1 by substituting the boojum texture and a circular 
boundary into Eq. (3.11). We have 



|V©o| 



iip fo( ek+iv 



) ~ e-^Ce*-*)] + [e^/^(e fc+ ^) + e^^e*"*)] } 



2 K / V+^)/ > fe -^) 
2na 2 

(1 - ae fc +^)(l - ac 



k — iip 



2na 2 f ™ k +*v 



1 — a 2 e 2k 



o e 



k—iip 



1 - ac k+l1 ? 1 
1 



1 



l 



^k—iip 



ae k-iip 
- 1 



+ 1 



(C.2) 



and 



a'( ( p-Q )Q ok + a"(>p-e )e 0lp 



ai r 
2 I 



-i(<p-e ) 



3 -i(¥>-e ) 



^k+iip jt^^k+iip^ 



-aie fe \e ie f~f» k+iv 



aiae 



7o(e fc 
1 

1 - ae k+l ^ 



)+c-^f' {e k -^)] 



■zk — itp 



(C3) 
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Substituting Eqs. (|CJ) and ( |C.3| ) into Eq. ( p.ll| ), all 

the if dependence cancels exactly. Equality is achieved 
by picking A = a\ae k — ao- We conclude that the circular 
domain with boojum texture is an equilibrium configu- 
ration for the case where p = 1. It is obvious that this is 
in general not true for any other p =/= I. 

We now examine the case of the bubble for p 1. We 

substitute the inverse boojum 0, into Eq. ( 3.11 ). We 
compute the following when c directors on the bound- 
ary favor pointing in the bulk, or pointing away from the 
center of the bubble, or a\ > 0, 



a'(n + <p- Qi)Q lk + a" {ti + tp - ©;)© 



(1 - ae- k+il ?) 2 (1 - ae~ k - i{ P) 2 



(C.4) 



When the c directors on the boundary favor pointing 
away from the bulk, or pointing toward the center of 
the bubble, or a\ < 0, we obtain 



a'(<p-Q i )e ik + a"( ( p-e i )e i¥ 



\ai\a 



- z tp 



(e- 



lip 



(1 



-k-\-ipy2 



(i 



-k—iip\2 



(C.5) 



And the following contribution is independent of the sign 
of III 



\ai\a 



1 



1 



-i — k—itp 



ae -k+i V 
l 



(C.6) 



Approximating 6 n w — & k , ©t « — V , and putting to- 
gether the above contributions, we find 



k" = - 



a a e k 

CLiOt 

a e k 

a e k 
2a\d 2 e 



1 



ae 



1 - ae~ k+i v (1 - ae- k+i f) 2 
e i2 ^ - 2ae~ k+ltp + 1 

c.c. 



c.c. 



(1 



ae" 



-k+iip\2 



(1 - a 2 e~ 2k y 2 * 



(1 



one 



2k 



ae -k+i V y 
e i2v 



c.c. 



(1 



ae" 



-k+iifi\2 



1 



(C.7) 



which reduces to Eq. (6.3). We keep in the expressions 
of k" only the apparently nontrivial terms and drop arbi- 
trarily the constant terms for convenience. In the actual 



evaluation of k, the constant is reinserted into k" to en- 
force the symmetry requirement k'(ip) = —k'(-ip). The 
boundary correction is taken to have little modification 
to the overall area enclosed and the treatment for fixing 
the area using the Lagrange multiplier A is ignored. 



APPENDIX D: EVALUATION OF THE 
THERMAL AVERAGES 



In this appendix, we describe the detail evaluation of 
(Q(cp)Q(ip')) o an d (exp[i^2qiQ(ipi)])o that we encoun- 
tered in Sec. |VII| . The average (O)o is taken with respect 
to 



h [Q] = \\ drivel 2 . 



(D.l) 



We first note that the extremum equation for © is given 
by Laplace's equation V 2 = in which the solution 
accommodates any boundary condition. Hence we can 
write Q(x,y) = @i(x,y) + ©2(2;, y) in general, where 



V 2 6i 



(D.2) 
(D.3) 



It can be shown further that ©i indeed minimizes ffo[0]. 
At low temperature, contributions from large ©2 are 
suppressed in the partition function. We shall assume 
©2 <C ©1 and i?o[@] is second order in ©2. The quan- 
tity ©2 can be neglected and we have, in the case of a 
domain, 



= 



OO 

E 

m— 1 



P m {a r , 



(D.4) 



Hq[0] can then be evaluated, 

00 

H = 4-7T V" m\A 



2 

m I ^ 



(D.5) 



where A m — a m Ro. 

We let &o(<p) = 0(i?ocos</?, i? siny>), 



(0o(^)0o(^)>o 
J(TO)e M9o(y')e-^ 
/(D9)e-^» 

Il„ / dA n dA* n J2 k J2k>( A ke lkv + Ale-^)(A k ,e* k 'v' + ,4*, e -« fc V) e - 2 ^E m A ™ A l 



k=l 



ik(<P-<p') 



2Y\ n fdA n dA* n e~ 2 ^ K ^ A ™ A " 
.„,„ V)1 f dA k dA* k A k Ale- 2 ^ kA «^ 
2 J dA k dA* k e- 2nf3KkAkA i 
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A ^-^ e ifc (¥>-¥>') _|_ Q-iH'-P-'-p') 



fe=i 
A 



— {in [l - e^-v') +ln l-e~ 4(v ~ v,) ]} 



■In 



4sin 2 --^-- 



-A In 



\x — x 



When ip — ► y/, we then introduce an ultraviolet cutoff i?o/a in the sum 

Ro/a 

(6 2 M) = A^I = Aln^. 



(D.6) 



(D.7) 



fc=i 



Evaluation of (exp[i <?i©(Vi)])o is best illustrated with 



/ e i[e( V i)-e( V2 )]\ o 



E E 

mi —0 m 2—0 



r ([ie(^)]^He^ 2 )] m2 ) . 



(D.8) 



We will apply Wick's theorem to compute the above average. We first note that the average vanishes when m\ + to 2 
is odd. When both m\ and mi are even, we have 

(lienor hom^o 

(e(^)e(^ 2 )>< (-e 2 (^))f^ (-e% 2 )) ^ 



mi!m 2 ! 



E 

i even 



(2m)!(2n 2 )! £ 

"o=0 



(-^-)!2 : 



2 1"" V 2 

2 "o / q2/,„ «o 



(9( yi )9( y2 )) 2 "° (-9 2 ( yi ))g 



-e% 2 )>r 



(2no)! 



(ni - n )!2 ni -™° (n 2 - n )!2"2- n o ' 



(D.9) 



where mi = 2ni, m 2 = 2n 2 , and / = 2uq. Similarly we can obtained for the cases where both m\ and m 2 are odd as 
below, 



(2ni + l)!(2n 2 + 1)! ^ 



n =0 



(e^pei^)) 2 "^ 1 (-e 2 (^- n " (-e 2 ( y2 ))^-» 

(2n + 1)! (ni - n )!2 n i-"o (n 2 - n )!2™2-"o 



(D.10) 



where toi = 2ri! + 1, m 2 = 2n 2 + 1, and I = 2n + 1. Combining these contributions in Eq. (D.8), we get 

/ e i[eOi)-e(<p s )]\ 

(e(yi)e( y2 ))g° (-e 2 ( yi ))^ (-e 2 (^ 2 ))^ 



E E E 

no=0 ni— n 2 =0 



«()! 



ni!2 r ' 



n 2 !2 f 



e (e( Vl )e( V2 ))o e 

[e( yi )-e(y 2 )] 2 



(e(yi)e(y 2 )>o (e( yi )e(y 2 ))o 



(D.ll) 



Although more tedious enumerations of the combinations 
of the correlation functions must be carried out in order 
to evaluate (exp[i qi&(ifi)])o, the same principle ap- 
plies. 

It can be observed in this simple c ase that the contri- 
bution in the partition function Eq. (7.2) of 



3 i[eOi)+e(v> a )h 



= c 



- 2 Aln-2-„-Aln- 



(D.12) 
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is approximately a factor c ln -. smaller than that of 
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